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' Abstract 

We present the results of an extensive lattice calculation of the renormaliza- 
tion constants of bilinear and four-quark operators for the non-perturbatively 0{a)- 
improved Wilson action. The results are obtained in the quenched approxima- 
tion at four values of the lattice coupling by using the non-perturbative RI/MOM 
renormalization method. Several sources of systematic uncertainties, including dis- 
cretization errors and final volume effects, are examined. The contribution of the 
Goldstone pole, which in some cases may affect the extrapolation of the renorma- 
lization constants to the chiral limit, is non-perturbatively subtracted. The scale 
independent renormalization constants of bilinear quark operators have been also 
computed by using the lattice chiral Ward identities approach and compared with 
those obtained with the RI-MOM method. For those renormalization constants the 
non-perturbative estimates of which have been already presented in the literature 
we find an agreement which is typically at the level of 1%. 



PACS numbers: 11. 15. Ha, 12.38.Gc, ll.lO.Gh. 



1 Introduction 



The increasing availability of computing resources has allowed a reduction of the statis- 
tical errors in current lattice QCD calculations at the level of few percent or less. The 
corresponding systematic uncertainties, however, are significantly larger, in some cases 
even by one order of magnitude. For this reason, a major effort in lattice calculations 
is aimed to reduce and to better control the various sources of systematic uncertainties. 
An important ingredient, in this respect, is the use of non-perturbative renormalization 
(NPR) techniques. By avoiding the perturbative expansion in the determination of the 
renormalization constants (RCs), these techniques remove the uncertainty associated with 
the truncation of the perturbative series and allow to reduce the systematic error involved 
in the renormalization procedure of the lattice matrix elements at a level of accuracy which 
is comparable or better than the present statistical one. 

The first NPR technique was proposed almost 20 years ago. Within this approach the 
RCs are determined by requiring that the renormalized lattice Green functions satisfy the 
proper chiral Ward identities (Wis) fij. This theoretically nice and numerically accurate 
approach suffers, however, of an important limitation: the use of chiral Wis only permits 
the determination of scale independent RCs (and mixing coefficients). After pp, it took 
about 10 years before other two NPR methods, which can be used in principle to compute 
any RC, were developed. They are based on the so called RI/MOM ^2] and Schrodinger 
Functional (SF) schemes [H]. In the RI/MOM scheme the renormalization conditions 
are imposed, non-perturbatively, on quark and gluon Green functions at large external 
momenta, while in the SF approach they are imposed on Green functions computed in a 
small volume. More recently, a new proposal for NPR has been suggested in |3], and it 
is based on the study of lattice correlation functions at short distance in x-space. The 
results of the first numerical investigations of this approach have been presented in ISHHl-^ 

The non-perturbative RI/MOM renormalization method has been already largely ap- 
plied in lattice calculations to determine the RCs of several classes of operators for different 
versions of the lattice action (see for example refs. [Hl-JSl)- In this paper we present the 
results of an extensive RI/MOM calculation of the RCs of the complete set of bilinear and 
four-fermion AF = 2 and AJ = 3/2 operators, for the non-perturbatively (9(a)-improved 
Wilson action J3] in the quenched approximation. The results are obtained at four values 
of the lattice coupling, in the range 6.0 < (3 < 6.45. Particular attention in the present 
calculation has been dedicated to the evaluation and control of systematic uncertainties: 

• Discretization errors have been evaluated by studying the behaviour of the RCs as a 
function of the lattice coupling, g{a), and of the renormalization scale. In the case of 
bilinear quark operators, since the action and the operators are non-perturbatively 
improved at C(a), leading discretization effects are of C(a^). These effects have been 
further reduced to O^g'^a'^) by using the results of a recent perturbative calculation 
of the relevant correlation functions jTHj. On the other hand, the four fermion 
operators considered in this study are not improved, and we are left in this case 
with leading discretization effects of 0{a). 

^See also ref. [7| for a related study in the context of the lattice non- linear ct- model. 
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• Finite volume effects have been examined by comparing the results of two indepen- 
dent simulations performed, at the same value of the lattice coupling (/3 = 6.0), on 
different lattice volumes. 

• Goldstone pole contributions: power suppressed contributions coming from the Gold- 
stone pole, which may affect the extrapolation to the chiral limit of the RCs of the 
pseudoscalar density and of the four-fermion operators coupled to the Goldstone 
boson, have been non-perturbatively subtracted. 

The scale independent RCs of the vector and axial-vector currents, Zy and Za, and the 
ratio Zp/Zs, have also been determined in this study by using the lattice chiral WI ap- 
proach, and the results are compared with those obtained with the RI/MOM method. 
We also compare our determinations of these constants with those obtained by the AL- 
PHA [T7j and LANL ^H] Collaborations by using the WI method, and our determi- 
nation of Zp with the one obtained by ALPHA [19j within the SF approach. We find 
an agreement which is typically at the level of 1%. Since the systematics involved in 
these approaches are different, these comparisons provide additional confidence on the 
high level of accuracy reached in the implementation of the RI/MOM non-perturbative 
method. 

Our final results for the RCs of bilinear quark operators and four-fermion operators 
are collected in tables I21IZI Preliminary results of the present study have been presented 
at the Lattice 2002 conference |H| I2()j. 

The plan of this paper is as follows. In sec. El we give the details of the numerical 
simulation and briefly review the non-perturbative RI/MOM method used to determine 
lattice RCs. The various sources of systematic errors involved in the calculation are 
discussed in sec. El where we also provide estimates of the corresponding uncertainties. 
The results for the RCs of bilinear quark operators obtained with the RI/MOM method 
are summarized in sec. El In this section, we also compare these results with those obtained 
by using the WI method, with predictions of one-loop boosted perturbation theory and 
with the results of refs. jni-CH]- The RI/MOM determination of the RCs of the four- 
fermion AF = 2 operators is discussed in sec. El and we end the paper by presenting 
our conclusions. More technical issues concerning the (9(a)-improvement of the RCs and 
the Goldstone pole contributions to the Green functions of four-fermion operators are 
collected in appendix A and B respectively. 

2 Details of the lattice calculation 

The lattice parameters used in this study are summarized in table ^ We have gener- 
ated 0(1000) gauge configurations in the quenched approximation, by using the non- 
perturbatively C(a)-improved Wilson action. The improvement coefficient csw has been 
determined in ref. ^3], as a function of the coupling constant. We have considered four 
different values of the lattice coupling, namely (3 = 6.0, 6.2, 6.4 and 6.45, corresponding 
to inverse lattice spacings varying approximately between 2 and 4 GeV. The estimates of 
the inverse lattice spacing given in table ^ and used in the numerical analysis, have been 
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/? = ^hl 




6.0 


6.0 


6.2 


6.4 


6.45 


csw 




1.769 


1.769 


1.614 


1.526 


1.509 


L^xT 




16^ X 52 


243 X 64 


243 X 64 


32^ X 70 


32^ X 70 


# conf. 




500 


340 


200 


150 


100 


a-\GeV) 




2.00(10) 


2.00(10) 


2.75(14) 


3.63(18) 


3.87(19) 


Ki 




0.1335 


0.13300 


0.1339 


0.1347 


0.1349 


K2 




0.1338 


0.13376 


0.1344 


0.1349 


0.1351 


«3 




0.1340 


0.13449 


0.1349 


0.1351 


0.1352 


/«4 




0.1342 




0.1352 


0.1353 


0.1353 






0.135225(5) 


0.135217(7) 


0.135815(3) 


0.135747(2) 


0.135686(2) 



Table 1: Summary of the lattice details and parameters used in this work. We also give the 
values of the inverse lattice spacing and of the critical hopping parameter (ncr)- 



determined from the study of the static quark anti-quark potential [21] by setting the 
scale ro = 0.530(25) fm, which corresponds to a~^(/3 = 6) = 2.0(1) GeV. The uncertainty 
takes into account the typical spread of the results obtained in the determination of the 
lattice spacing in the quenched theory. We emphasize that mainly ratios of scales, rather 
than their absolute values, are relevant for the present calculation and for these ratios the 
static quark anti-quark potential provides an accurate determination. 

For each value of the coupling, the RCs have been computed at four values of the 
light-quark mass and eventually extrapolated to the chiral limit. ^ The values of the 
Wilson hopping parameter used in the simulations are listed in table ^ together with 
the corresponding critical values determined from the vanishing of the axial WI mass. 
The Wilson parameters correspond to bare quark masses which, in lattice units, are 
in the range 0.01 < aruq < 0.05. The same set of gauge configurations and quark 
propagators used in this study have been also used in ref. |^, where additional details on 
the simulation can be found. In order to study finite volume effects we have considered two 
independent lattice simulations at /3 = 6.0, performed on different lattice sizes. The main 
purpose of the simulation on the larger volume was the study oi K —>■ ttit decays fl^ |23 
and this motivates the different choice of the values of quark masses. 

The non-perturbative determination of the RCs with the RI/MOM method is based 
on the numerical evaluation, in momentum space, of correlation functions of the relevant 

^The RI/MOM is, by definition, a mass independent renormalization scheme. 
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operators between external quark and/or gluon states. In the case of the bihnear quark 
operators Or = qTq, where T = S, P,V, A,T stands respectively for /, 75, 7^, 7/^75, a^i^, 
the relevant Green function is 

Gr{p,p') = f d'xd'ye-'P-'+'P' y (g(x)Or(0)^(i/)) , (1) 



where q and q are renormalized quark fields. By using Gr{p,p') and the renormalized 
quark propagator, 

S{p)= /rf^xe-*P-^(g(x)f(0)), (2) 



one evaluates the forward amputated Green function 

AT{p) = S{pr'Gr{p,p)S{p)-'. (3) 

The RI/MOM renormalization method consists in imposing that the forward amputated 
Green function, computed in the chiral limit in a fixed gauge and at a given (large) scale 
p'^ = /x^, is equal to its tree-level value. In this study we work in the Landau gauge, but 
different choices of generic covariant gauges have been also considered in the literature . 
In practice, the renormalization condition is implemented by requiring 

Zr Trip) = Zr Tr [Ar(p) Pr] |p2=^2 = 1 (4) 

in the chiral limit, where Pr is a Dirac projector satisfying Tr [FPr] = 1. The RC of the 
quark field, which enters the definition of the renormalized quark propagator, is obtained 
by imposing the condition 



^Sip)-' 



p2 



Za — Tr 
^ 12 



i>sip)-' 



p2 



(5) 



in the chiral limit. In the numerical calculation, throughout this study, we adopt the 
definitions 

1 ^ 1 ^ 

= — V sin^(ap^) , / = - V 7,.sin(ap^) (6) 

and, in order to minimize discretization effects, we select momenta with components p„ = 
(27r/L^) in the intervals = ([0, 2], [0, 2], [0, 2], [0, 3]) and = ([2, 3], [2, 3], [2, 3], [4, 7]) 
(Ljy is the lattice size in the direction u). 

The renormalization scale fi, introduced in eqs. (j^ and must be much larger 
than Aqcd, in order to be able to connect the RI/MOM renormalization scheme to other 
schemes by using perturbation theory, and much smaller than the inverse lattice spacing, 
to avoid large discretization errors: 

Aqcd < /i < 7r/a . (7) 

An important point to be mentioned is that the RCs of bilinear quark operators, ob- 
tained with the RI/MOM technique in the chiral limit and at sufficiently large values of 
the external momentum, are automatically improved at 0{a). The reason is related to 
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M-2 



M-2 




Figure 1: The scale independent RCs Zy , Za and the ratio Zp/Zs (left) and the RGI combi- 
nations Zp'^^ , Zg*^^ and Z^'-'^ (right) are shown as a function of the renormalization scale, at 
(3 = 6.2. The solid lines represent the results of a constant fit to the data. 



chiral symmetry and to the fact that 0{a) on-shell improvement corrections, proportional 
to the coefficients cy, ca and ct, vanish for Green functions computed between forward 
external states, at zero momentum transfer. An explicit proof of this statement is given 
in appendix A. Being the improvement coefficients known with some uncertainty, the pos- 
sibility of neglecting them in the calculation improves the accuracy of the determination 
of these RCs. In contrast, on-shell (9(a)-improvement of the operators is required when 
the RCs are computed by using the WI approach. In this case, in order to improve the 
vector and axial-vector current operators, we use the values of the coefficients cy and ca 
determined non-perturbatively in refs. QB] ■ 



3 Systematic errors 

In this section we discuss the RI/MOM non-perturbative evaluation of the RCs of bilinear 
quark operators by examining in details the various sources of systematic errors and 
providing estimates of the related uncertainties. 

3.1 Renormalization scale dependence and discretization effects 

The results for the scale independent RCs Zy and Za, and for the ratio Zp/Zs, obtained 
with the RI/MOM method are presented in the left plot of fig. ^ as a function of the 
renormalization scale. We show, as an example, the results obtained aX (3 = 6.2. The 
good quality of the plateau indicates that 0(a^) discretization effects are rather small, at 
the level of the statistical errors, even in the region a/x ^ 1 considered here. 
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For the scale dependent bilinear operators, O = S,P, T, we show in the right plot of 
fig. H the renormalization group invariant (RGI) combinations 



Z^^^=Zo(/i)/Co(^), (8) 

where the evolution function Cc)(/i) = exp J"'^^^ da'jo{o:)/P{o!) , with P{a) and 7c>(«) 
the beta function and the anomalous dimension of the relevant operator, is introduced 
in order to explicitly cancel, at a given order in perturbation theory, the scale depen- 
dence of the RCs. In the RI/MOM scheme, these functions are known at the N^LO for 
Zt [2ni and at the N^LO for Zs and Zp j2ZI. In their numerical evaluation, we use the 
determination of the strong coupling constant obtained, in the quenched theory, from 
A^~°'' = 225(20) MeV ^Hl- From the quality of the plateau shown in the right plot of 
fig. In discretization effects appear to be very small, even at (unexpectedly) large values 
of a/i.^ 



3.2 Perturbative correction of O{g^o?) discretization effects 

Any deviation from the predicted constant behaviour at large momenta, as observed in 
fig- E and in analogous results obtained at other values of (3, signals the presence of either 
0{a?) discretization effects or of higher-order perturbative corrections not included in the 
evaluation of the evolution function Co{n)- These deviations are found to be larger on 
the coarser lattice, corresponding in our case to (3 = 6.0. One can also notice in fig. Q 
that the quality of the plateau is worse in the case of the RC of the scalar density Zs- As 
discussed below, we interpret this to be due to the presence of larger discretization errors. 

The leading 0{g'^a?) discretization errors in the lattice expressions of the amputated 
Green functions Ar(p) of eq. Q have been recently computed by using lattice perturba- 
tion theory [IS^. In order to further investigate the effect of discretization errors in the 
calculation of the RCs, we use the results of [TH] and subtract the perturbative 0{g'^a^) 
contributions from the amputated Green functions Ar(p) computed non-perturbatively in 
the numerical simulation. 

The most sensitive to the subtraction of the 0{g'^a^) terms is the constant Zs- The 
correction decreases the final estimate of this RC by approximately 5% and we find that 
after the subtraction the quality of the plateau in the RGI combination is significantly 
improved. The result is shown in the left plot of fig. |2l at /? = 6.2. On the other hand, 
the effect of the 0{g'^a?) correction is found to be negligible for Zp, as shown in the right 
panel of fig. El In the case of the RCs Zt the correction amounts to approximately 3% 
and for Zy and Za to approximately 1%. 

All the results for the RCs of bilinear quark operators discussed in the following will 
be those with the terms of O{g^o^) subtracted away. After this correction, the systematic 
error due to residual 0{a?) effects is evaluated from the quality of the plateau of the RGI 
combinations as a function of the renormalization scale. This error will be included in 
the final evaluation of the RCs presented in sec. 0] 

■^The Goldstone pole contributions to the ratio ZpjZg and to the combination Zp"-^^ shown in fig-H 
have been non-perturbatively subtracted as discussed in sec. 13.41 
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Figure 2: The RGI combinations Zg'^^ (left) and Zp^^ (right) as obtained with and without 
the subtraction of the 0(0^0?) contributions computed in perturbation theory. 

3.3 Comparison of results at different values of the coupling 

Another possibility to study discretization effects is based on tlie comparison of tfie RCs 
computed at different values of the lattice spacing. Though the RCs are obviously func- 
tions of the lattice spacing (the UV cutoff in the lattice regularization) , their dependence 
on the renormalization scale should be universal (i.e. independent of a), and only fixed 
by the anomalous dimension of the corresponding operators. Therefore, the dependence 
on the renormalization scale should cancel in the ratio 

R{a,a) = Z{a, fi)/Z{a, fi) , (9) 

up to discretization effects. In fig. El we show as an example the results for the RCs 
Zy and Zp obtained, at the four values of (3, as a function of the renormalization scale. 
Each RC, Z{a,fj,), has been rescaled by the factor R{a,a) of eq. ©, where we have 
chosen as a reference scale a the value of the lattice spacing at /? = 6.4. From fig. |2l we 
observe that the results obtained at the different values of /3 all lies, for large values of the 
renormalization scale, on the same universal curve, thus confirming that discretization 
effects are well under control within the statistical errors. The figure also shows that 
the renormalization scale dependence of the RCs is in very good agreement with the one 
predicted by the anomalous dimensions of the relevant operators. On the other hand, we 
note that at small values of the renormalization scale the results obtained at different (3s 
show some disagreement, and that the dependence on the renormalization scale differs 
from the one predicted in perturbation theory. These discrepancies may be due to finite 
volume effects, higher order perturbative corrections and also to power suppressed non- 
perturbative contributions which might not be negligible in the region of small momenta. 
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Figure 3: Values of Zy (left) and Zp (right) as obtained at the four values of the lattice coupling 
as a function of the renormalization scale. The results have been rescaled by the factor R{a, a) 
defined in eq. (0). The solid lines represent the scale dependence predicted by the relevant 
anomalous dimensions (7 = in the case of Zy ). 

For this reason, the results obtained at small momenta are not considered in our final 
determination of the RCs. 



3.4 Subtraction of the Goldstone pole 

The validity of the RI/MOM approach relies on the fact that non-perturbative contribu- 
tions to the Green functions vanish asymptotically at large p^. Among the bilinear quark 
operators, however, specific care must be taken in the study of the pseudoscalar Green 
function Vp since in this case, due to the coupling with the Goldstone boson, the leading 
power suppressed contribution is divergent in the chiral limit j21l2Hll2ni- The operator 
product expansion (OPE) of Fp reads in fact: 

Fp(p2, m) ^ ci{p\ m) + C2(p^ m) ^ + , (10) 

where m is the quark mass, (gg) is the quark condensate in the chiral limit and ci^2 
are Wilson coefficients. In practice, since we work in a region of light quark masses 
(m^/p^ -C 1) varying in a short interval, the mass dependence of the Wilson coefficients 
can be neglected. We extract the RC Zp from the short- distance behaviour of Fp at large 
p^, represented in eq. ()1U|) by the coefficient Ci. In order to determine the contribution of 
the Wilson coefficient Ci we have then followed two different approaches. 

In the first one the Goldstone pole contribution proportional to C2 has been subtracted 
by constructing the combinations [SU] : 

rSUB/ 2 N miFp(p2 mi) -m2Fp(p2 ma) 

F|,^ (p,mi,m2) = , (11) 

mi — m2 
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Figure 4: The RGI combination Zp as obtained with and without the non-perturbative 
subtraction of the Goldstone pole, at j3 = 6.4. 



where mi and m2 are non degenerate quark masses.^ The effect of the subtraction on 
the RGI combination Zp*^^ is shown in fig. 0] Though both the unsubtracted and the 
subtracted determinations of Zp^^ converge to approximately the same value at large p^, 
a clear plateau is never observed in the first case, at least not in the region of momenta 
explored in this study. For this reason, we have derived our first determination of the RC 
Zp by using the subtracted Green function FfP^ defined in eq. (fTTj) . 

In the second approach we first study the dependence of Fp on the quark mass |2H! • 
Guided by the OPE of eq. (jlUj) . we have performed a fit of Fp to the form^ 

Fp(/,m) = + + (12) 

m 

The singular Goldstone pole contribution to Fp is represented in this expression by the 
term proportional to B{p'^). The coefficient ^(p^) provides instead, at large p^, the short- 
distance asymptotic behaviour of the Green function in the chiral limit. From this coef- 
ficient we have therefore obtained our second determination of Zp. 

We find that the RC determined with the two approaches, namely either from the 
subtracted Green function Fp^^ in eq. or from the coefficient A(p^) of eq. ()12|) . are 
in perfect agreement, the difference being less than 1% at all values of the coupling (3. 

In order to further investigate the Goldstone pole contribution to the pseudoscalar 
Green function, we have studied the momentum dependence of the coefficient -B(p^) by 

''in this analysis we used the vector WI definition of the quark mass, iriq = ^ (I/k, — I/kct)- We also 
verified that using the axial WI definition leads to practically indistinguishable results. 

^ Since the vertex Green function has been computed by inserting operators with both degenerate and 
non-degenerate quark masses, the 3-parameter fit in eq. (|12() is performed on a set of 10 data points. 
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performing a fit to the form 

5(p')=« + 4 + ^- (13) 

We find that in the region of relatively large external momenta (p > 2 GeV) this fit 
provides a good description of the lattice data. The values of the coefficient a, which 
may be only generated by pure lattice artifacts, are found to be consistent with zero at 
all values of the lattice coupling. From the results for the coefficient (3 we have derived 
an estimate of the quark condensate, which turns out to be in agreement with other 
lattice determinations of the same quantity based on different approaches [Sj-ISEl- The 
details of the extraction of the quark condensate from the pseudoscalar Green function 
are presented in a separate publication ^Tj . 



3.5 Chiral extrapolations 

The renormalization condition in eq. (^, which defines the RI/MOM renormalization 
scheme, must be implemented in the chiral limit. This ensures that RI/MOM is a mass 
independent renormalization scheme. Since the non-perturbative determinations of the 
RCs are obtained in the calculation at non vanishing values of the quark masses, a final 
extrapolation of the results to the chiral limit must be eventually performed. 

After the contribution of the Goldstone pole has been non-perturbatively subtracted 
from the pseudoscalar correlator, all Green functions of bilinear quark operators are ex- 
pected to be smooth functions of the quark masses. In the region of masses considered in 
this study, this mass dependence is found to be rather weak and consistent with a simple 
linear behaviour. For this reason, the extrapolation to the chiral limit is easily performed. 
For illustrative purposes, we show in fig. El the chiral extrapolation of the RCs Zy and Zs 
at the four values of the lattice coupling. The smooth dependence on the quark masses 
suggests that the systematic error introduced by the chiral extrapolation is smaller than 
the statistical error and can be therefore safely neglected. 



3.6 Finite volume effects 

In order to study finite size effects, two independent simulations on different lattice sizes 
at /3 = 6.0 have been considered. The smallest size corresponds approximately to the 
same physical volume used at the other values of f3. From this analysis we find that the 
differences between the values of the RCs obtained at the same value of the renormalization 
scales from lattices of different sizes are smaller than the statistical errors. For this reason, 
we have not included an additional systematic error due to final volume effects in our 
final determination of the RCs at all values of the lattice coupling. An example of this 
comparison, for the RCs Zy and Zs, is shown in fig. El 



4 Results for the RCs of bilinear quark operators 

In this section we present our final results for the RCs of bilinear quark operators as 
obtained with the RI/MOM method. We then compare these results with those obtained 
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Figure 5: Values of Zy (left) and Zs (right) as obtained at the four values of the lattice 
coupling as a function of the bare quark masses. The solid lines represent the results of a linear 
extrapolation to the chiral limit. 
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Figure 6: Comparison of Zy (left) and Zs (right) at (3 = 6.0 as obtained on lattices of different 
sizes as a function of the renormalization scale. 
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RI/MOM 


WARD IDENTITIES 


BPT 1-loop 






This work 


This work 


LANL 


ALPHA 


csw = 1 






6.0 
6.2 
6.4 

6.45 


0.772(2)(2) 
0.783(3) 
0.801(2) 
0.803(3) 


0.774(4) 
0.789(2) 
0.804(2) 


0.770(1) 
0.7874(4) 
0.8018(5) 


0.7809(6) 
0.792(1) 
0.803(1) 
[0.808(1)] 


0.7820 
0.7959 
0.8076 
0.8103 


0.8504 
0.8463 
0.8480 
0.8488 


Za 


6.0 
6.2 
6.4 
6.45 


0.812(2) 
0.819(3) 
0.832(3) 
0.833(3) 


0.856(17)(15) 
0.812(5)(2) 
0.843(10)(1) 


0.807(8) 
0.818(5) 
0.827(4) 


0.791(9) 
0.807(8) 
0.827(8) 
[0.825(8)] 


0.8038 
0.8163 
0.8269 
0.8293 


0.8693 
0.8624 
0.8628 
0.8633 


Zp 
Zs 


6.0 

6.2 


O.o/0(4j(oj 

0.877(5) 


0.893(20)(19) 
0.877(5)(1) 


0.842(5) 
0.884(3) 


[0.840(8)] 
[0.886(9)] 


0.9398 
0.9449 


0.9444 
0.9545 




6.4 
6.45 


0.894(3) 
0.897(4) 


0.914(10)(1) 


0.901(5) 


[0.908(9)] 
[0.912(9)] 


0.9491 
0.9500 


0.9594 
0.9603 



Table 2: Values of the scale independent RCs as obtained with the RI/MOM method, the WI 
method and one loop boosted perturbation theory (BPT). 

by using the WI approach (for the scale independent constants), one-loop boosted per- 
turbation theory and with other results presented in the literature. 

4.1 RI/MOM results 

Our final estimates of the RCs of bilinear quark operators obtained with the RI/MOM 
method are derived by fitting to a constant the RGI combinations shown in fig. ^ (and 
similarly at the other values of (3) after having corrected the 0{g'^a?) discretization errors, 
as discussed in sec. 13.21 and performed the extrapolation to the chiral limit. The results for 
the scale independent and scale dependent RCs are collected in the third column of tablesEl 
and El respectively. The second error, when quoted, represents the systematic uncertainty 
due to residual discretization effects, estimated from the quality of the plateaus of the 
RGI combinations. 
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RI/MOM 


SF 


BPT 1-loop 






This work 


ALPHA 


Csw — 1 




Zf(l/a) 


o.U 
6.2 
6.4 
6.45 


n Q'iCi( '~)\( f\\ 
u.ooy^^z )\P) 

0.850(2) 

0.861(2) 

0.860(3) 





0.8151 
0.8269 
0.8369 
0.8391 


0.8821 
0.8751 
0.8750 
0.8754 


Zf(l/a) 


p. n 

D.U 

6.2 
6.4 
6.45 


0.642(3) 
0.671(4) 
0.680(7) 




0.6685 
0.6896 
0.7074 
0.7115 


0.6257 
0.6558 
0.6794 
0.6846 


Z#^(l/a) 


6 n 

6.2 
6.4 
6.45 


0.564(4) 
0.600(4) 
0.610(8) 


[0.57(1)] 
[0.60(1)] 
[0.61(1)] 


0.6248 
0.6487 
0.6689 
0.6735 


0.5877 
0.6236 
0.6498 
0.6555 


Z«^(l/a) 


6.0 
6.2 
6.4 
6.45 


0.881(2)(3) 
0.876(2) 
0.884(3) 
0.883(2) 




0.8417 
0.8518 
0.8603 
0.8622 


0.9442 
0.9259 
0.9190 
0.9180 



Table 3: VaJues of the scale dependent RCs as obtained with the RI/MOM method, the SF 
method and one loop boosted perturbation theory (BPT). The results are expressed in the 
RI/MOM scheme at the scale /x = 1/a. 
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4.2 Comparison with results from Wis 

In order to compare the results obtained with the RI/MOM method with those obtained 
by using a different non-perturbative approach, we have also computed the scale indepen- 
dent RCs Zy, Za and the ratio Zp/Zs by studying the lattice chiral Wis [HEH]- We have 
performed the calculation at three of the four values of the lattice coupling considered for 
the RI/MOM study, namely (3 = 6.0, 6.2 and 6.4. The vector and axial vector current 
operators have been improved at 0{a) by using the values of the coefficients cy and ca 
determined in refs. |14| I18j. 

The RC Zy of the local vector current has been determined by imposing 

2p E E {P{x)Vo{y)pm = ^ E {My)Pm , (m) 

- y y 
where 2p = (E^ Vo(Ao(2;)P(0)))/(Ea?(^(a;)P(0))). We find that an independent deter- 
mination based on the WI 

2pEE (^(^)V^)(l/)^o(0)) = ^ E (^o(y)Ao(O)) , (15) 

- y y 

provides results consistent with those obtained from eq. (fT^ but with larger statistical 
errors. 

The RC of the axial current Za and the ratio Zp/Zs have been determined by using 
the fiavour non-singlet identities 

2pE E {P{x)VMAm) = -%Y. ^yAy)ym + ^ E ^My)Mm , (16) 

X y Ay V y 

with k summed over 1,2,3, and 

2pEE (^(^)^(2/)^(o)) = E {Piy)Pm + E {S{y)sm (17) 

X y y y 

respectively. 

The results are presented in the fourth column of table El The first error quoted in the 
table is statistical while the second one represents the systematic uncertainty introduced 
by different estimates of the value of the improvement coefficient ca- At /? = 6.0, in 
particular, the determination ca = —0.037 obtained in ref. |TH] is in disagreement with 
Ca = —0.083 estimated in ref. J3] by approximately a factor two. This difference in- 
troduces in turn a large uncertainty in the evaluation of Za and even more of the ratio 
Zp/Zs, as can be seen from the systematic errors quoted in tabled At the larger values 
of (3 the estimates of ca obtained in refs. ^3] and JHl are in better agreement and the cor- 
responding uncertainty in the evaluation of the RCs with the WI method is significantly 
reduced. 

Our final results for the RCs of bilinear quark operators, as obtained from the RI/MOM 
and the WI methods (the latter for the scale independent constants) are shown in fig. [3 
where they are also compared with the predictions of 1-loop boosted perturbation theory. 
We observe a very good agreement, typically within one standard deviation, between 
the RI/MOM and the WI determinations. The central values obtained with the two 
approaches differ by less than 1% for Zy, 5% for Za and 3% for Zp/Zs- 
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Figure 7: Values of the RCs as obtained from the RI/MOM method (fiUed circles), the WI 
method (empty circles) and 1-loop boosted perturbation theory with csw = 1 (solid lines). 



4.3 Comparison with perturbation theory 

The expressions of the RCs of bihnear quark operators computed in perturbation theory 
(PT) [SH] are functions of both the couphng constant g'^ and of csw, the improvement 
coefficient of the fermionic action. In the numerical calculation csw has been fixed to 
its non-perturbative value ^3]. In PT one can use either this non-perturbative estimate 
or the perturbative expression of csw, consistently evaluated at the proper order in g"^. 
The difference between the two evaluations represents an uncertainty which is of the same 
order of the other higher order terms neglected in the perturbative expansion. 

In tables |21 and El we present the predictions for the RCs obtained by using one-loop 
boosted PT in two cases, with csw fixed either to its tree-level value ^, csw = or to its 
non-perturbative estimate. The differences between the two sets of results, in the range 
6.0 < 13 < 6.45, are roughly of the order of 10%. This can be considered therefore the 
typical size of uncertainty associated with the predictions of one-loop boosted PT. The 
perturbative estimates, shown for illustration in fig. [71 are those obtained with csw = 1- 

The comparison between the non-perturbative RI/MOM determinations and the cor- 
responding perturbative estimates at one loop, presented in tables El and El and illustrated 
in fig. shows that the differences are indeed at the expected level of 10%. A somewhat 
larger discrepancy is observed for Zp. In this case, at /? = 6.0, it ranges between 10% 
and 20%. As expected, we find that in all cases the agreement between perturbative and 
non-perturbative determinations of the RCs improves by going towards smaller values of 



^In the perturbative expansion of the RCs csw only enters at 0{g^). 
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the coupling g^. 

We conclude this discussion by observing that the differences between perturbative 
and non-perturbative results for the RCs are significantly larger than the typical size 
of statistical errors in present lattice calculations. This is the reason why the use of 
NPR techniques should be considered, at present, a fundamental ingredient in the lattice 
determinations of the physical hadronic matrix elements. 

4.4 Comparison with other results in the Uterature 

The scale independent RCs Zy^ Za and the ratio Zp/Zs, for the C(a)-improved Wilson 
action considered in this paper, have been also computed by the ALPHA [T7j and 
LANL Collaborations by using the WI method. The results are collected in table |21 
We have quoted in square brackets the estimates of the RCs which are not given directly 
in the original papers but can be inferred from the published results. Specifically, the 
ALPHA results for the RCs at the value of the coupling (3 = 6.45 have been obtained 
by using the fits of these quantities, in terms of rational functions of g"^, performed in 
refs. |16l|T7j. We have also quoted in square brackets in table E] the ALPHA results for 
the ratio Zp/Zs, which have been obtained by combining the results for the quantity 
Zp/ {Za Zs) presented in ref. jTJj with the estimates of Za given in ref. [TT)] . 

The comparison between the RI/MOM results presented in this paper with those 
obtained by ALPHA and LANL by using the WI method shows an agreement which is, 
in most of the cases, at the level of 1% or even better. Slightly larger differences, at the 
level of 3%, are only observed for the estimates of Za and Zp/Zs at = 6.0. 

In tableiniwe also present (in square brackets) the results for the RC of the pseudoscalar 
density Zp obtained by the ALPHA Collaboration ^Hj by using the SF approach. Contin- 
uum perturbation theory at one loop has been used in this case to convert these estimates 
from the SF to the RI/MOM renormalization scheme, at the proper value of the renormal- 
ization scale. Also in this case, we find an agreement which is at the level of 1% or better. 
In conclusion, we believe that these comparisons, performed among results obtained by 
different groups and by using different approaches, provide strong evidence of the high 
level of accuracy reached at present by the several NPR methods and, in particular, by 
the RI/MOM approach. This also provides us with additional confidence on the quality of 
the results obtained for those RCs, like Zs and Zt, for which non-perturbative estimates 
have not been presented before with methods different from RI/MOM. 

5 Four-fermion operators 

We now come to the case of AF = 2 and A/ = 3/2 four-fermion operators which, 
in a mass independent renormalization scheme, share the same set of RCs. From a 
phenomenological point of view these operators play an important role since they enter 
the effective Hamiltonian of weak interactions. Matrix elements of AF = 2 operators, 
for instance, control the — K'^ and — mixing amplitudes, in both the Standard 
Model and its SUSY extensions (see for instance [101 E])- Operators with A J = 3/2 
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provide the amplitude of K ^ tctt decays in the 1 = 2 channel, and enter therefore the 
theoretical estimates of the A/ = 1/2 rule and of the direct CP violation parameter e'/e. 

Due to the explicit breaking of chiral symmetry induced by the Wilson term, four- 
fermion AF = 2 and AI = 3/2 operators mix with operators of the same dimension but of 
different "naive" chirality PJ 02] • ^ In order to define the complete basis of dimension-six 
four-fermion operators which mix under renormalization we introduce the notation 



1 



[(9irg2)(g3r'94) ± (girg4)(g3r'92)] • 



and 0[rir2±r3r4] = 0[rir2] ^^[r^r^]- The basis consists then of ten parity conserving (PC) 
and ten parity violating (PV) operators 



Qf 


— '^[VV+AA] 


Qf 


— ^[VA+AV] 


Qf 


— '^[VV-AA] 


Qf 


— ^[VA-AV] 


Qf 


— '^[SS-PP] 


Qf 


— '^[SP-PS] 


Qf 


— '^[SS+PP] 


Qf 


— ^[SP+PS] 


Qf 


= 


Qf 


= 0±-. 

TT 



(19) 



As extensively discussed in |43^, with Wilson fermions the renormalization pattern of 
the PV sector follows the "continuum" one 
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(20) 



In the PC sector, instead, the explicit breaking of chiral symmetry induces an additional 
mixing which is parameterized in terms of a matrix A, 
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In matrix notation, we can write 
where / is the 5x5 unit matrix. 



± 



Q^ = Z±[J + A±]Q 



± 



(21) 



:22) 



''In the SU{2) isospin limit there is no mixing with operators of lower dimension. 
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5.1 The RI/MOM method for four-fermion operators 

The procedure implemented to compute non-perturbatively the RCs of four-fermion op- 
erators within the RI/MOM scheme is a generahzation of that explained for bilinears. It 
has been presented in ref. [43]. Here we briefly outline this procedure and refer the reader 
to that reference for more details. The calculation proceeds in a number of steps: 

1. For each four-fermion operator Oi, we start by computing the four-point Green 
function in the Landau gauge 

Gi(xi, X2, X3, X4) = ( gi(xi) %{x2) Oi{0) qsix^) "q^ix^) ) , (23) 

where q and q are renormalized quark fields. For convenience, we take the Fourier 
transform of Gi with all external quark legs at equal momentum p. In this way 
we obtain Green functions in momentum space of the form Gi{p)'^^!^^, where the 
superscripts and subscripts are respectively the colour and spinor indices of the four 
external fields in eq. (jSBI)- 

2. Green functions are then amputated by multiplying Gi{p) by four renormalized 
inverse quark propagators S'^{p): 

HpT^'^^'s' = S-\p)± S-\prX G^ipfXs S-\P)%' S-\p)f; . (24) 

3. Projection operators P^, with the label i running over all the operators of the basis, 
are introduced which satisfy the orthogonality relation 

Tr{A?^°^P,^} = 5,u. (25) 

In the previous equation, the trace is taken over colour and spinor labels and A^'^^'' 
is the tree- level amputated Green function of the operator Oi. The explicit form of 
the projection operators is given in ref. [43] . 

4. Finally, the RI/MOM renormalization conditions are implemented. In the case of 
PV operators these conditions take the form 

^5 rj(p)| . = Zt, Tr {^f{p)P^}L_. = , (26) 

while for the PC ones they read 



(J + A±),, T%ip) ^ Zt (/ + A±),, Tr {Af (p)P,^} = 5., . (27) 

In order to ensure the mass independence of the renormalization scheme, the ma- 
trices Z and A are determined after the extrapolation to the chiral limit of the 
correlation functions entering eqs. and (j^Zj). 
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5.2 Subtraction of the Goldstone pole 

As discussed in sec. 13.41 for the case of the bihnear pseudoscalar operator, a possible 
difficulty in the implementation of the above procedure comes from the coupling of the 
operators to the Goldstone boson. In the renormalization of the four-fermion operators 
this difficulty was already observed in ref. pH] . 

By using the LSZ reduction formula we show in appendix B that in the PV case 
only a single pion pole can be present, while in the PC case the appearance of a double 
pole is also possible. This is the main difference with respect to the case of the bilinear 
pseudoscalar operator. 

In the PV case, by performing the OPE we find that the dependence on the quark 
mass m and the external momentum p of the projected Green function F can be written 



as^ 



where dots represent non-perturbative contributions which vanish at large momenta as 
or faster. In eq. (j^Hj) . T is the short-distance form factor, from which we want to 
extract the matrix Z of RCs, while the term proportional to AF is the contribution from 
the Goldstone boson propagator. An example of this behaviour is shown in the left plot 
of fig. IHl where the Green function F33(m,p^) (the operator Q3 happens to be strongly 
coupled to the Goldstone boson) is plotted as a function of 1/m for different values of p^. 
A linear dependence of F33 on 1/m is clearly visible in the plot, with a slope which is a 
decreasing function of the external momentum. 

In order to subtract the Goldstone boson contribution and to compute the RCs directly 
from the value of F in the chiral limit we follow the same strategy discussed in sec. 13.41 for 
the bilinear pseudoscalar operator. This strategy is based on two different approaches. 

In the first case we study the dependence on the quark mass of the projected Green 
functions Fjj by performing a fit to the form 

F,,(m,p2) = A,,{p') + + m . (29) 

m 

The coefficient A{p^) gives a determination of F in the chiral limit, whereas B{p^) is asso- 
ciated with the Goldstone boson contribution proportional to AF. The term proportional 
to C{p'^) comes from the expansion of F and AF in the quark mass. The results of the fit 
to eq. (j2ni); in the case of the correlation function F33, are shown as solid lines in the left 
plot of fig. IHl The matrix Z of RCs is then extracted from the coefficient A{p'^) at large 
p^. 

As an additional check of this analysis, we study the dependence of the coefficient 
B{p^) on the external momentum. In the right panel of fig. IHlwe plot i?33(p^) at /5 = 6.2 
as a function of l/{apY for {apY > 0.6. The solid line in the figure represents the result 
of the fit to the form 



*For notational simplicity, we drop from now on, unless necessary, the plus or minus superscript. 
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r22(m, p ) at P = 6.2 
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Figure 8: Left: The projected Green function T^^^mjp^) as a function ofl/am for five different 
values ofp^. The sohd hnes represent a fit to eq. Right: the coefRcient Bs3{p'^) of eq. |23) 

as a function of l/{apY and the corresponding fit to eq. if^) . 



We find that, within the errors, eq. (j30|) gives a good description of the data with a = 
—0.0009(44), a result which is well compatible with zero. This is consistent with the 
expectation that the Goldstone pole can only appear in power suppressed non-perturbative 
contributions. 

The second approach used to subtract the Goldstone pole is the one suggested in 
ref. ISni and also discussed in sec. 13.41 for the bilinear pseudoscalar operator. Applied to 
the case of the four-fermion operators, this method consists in eliminating the contribution 
of the Goldstone boson by constructing the combinations 



mi Tij{mi,p'^) -1712 Tij{m2,p'^ 
mi — m2 



(31) 



The subtracted correlation functions are then extrapolated to the chiral limit by perform- 
ing a fit to the form 



rrK,^2,/) 



Aijip"^) + Cij{p^){mi + m2) . 



(32) 



We find that the results for A(j9^), obtained from a fit which uses only the four lightest 
values of quark masses, differ from those of A^p"^) of eq. (j29|) by less than 1% in all 
the range of considered in this study. Therefore, at this level of accuracy, the two 
approaches provide the same estimate of the RCs. 

In the PC case also a double Goldstone boson pole can be present. For this reason, 
we fit in this case the subtracted correlation function of eq. ()31|) to the form 



rf7^(mi, m2,p^) = Aj{p^) + Cij{p^){mi + ma) + 



mim2 



(33) 
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2 2 

|j, a 

Figure 9: A34 as function of fi at P = 6.2 in three cases: without the subtraction of the Goldstone 
pole (circles), with the subtraction of the single pole (squares) and with the subtraction of both 
the single and the double poles (diamonds). 



which differs from eq. ()32j) for the presence of the last term. As before, we compute 
the RCs from the coefficient Had we used this form in the PV case, a value of 

E{p^) fully compatible with zero would have been obtained. In the PC case instead, our 
numerical results clearly show the presence of such contribution, which is however much 
smaller than the single pole contribution proportional to -B(p^). 

The effect of single and double poles, when not subtracted, is clearly visible in several 
matrix elements, both of the scale dependent RCs, Z(^a,gQ) and Z(fia,gQ), and of the 
scale independent mixing coefficients A((7q). An example of the latter case, namely the 
one of the matrix element A34, is shown in fig. El It is very reassuring that A34 becomes 
almost independent of the scale, as should be the case, once the pole contribution has 
been eliminated. For some other matrix elements, and in particular in the case of Zu and 
of the corresponding mixing coefficients Ai^ which are relevant for the lattice estimates 
of — and — mixing in the Standard Model, the effect of the subtraction of 
the Goldstone pole is absolutely negligible. 

5.3 Renormalization scale dependence, discretization and finite 
volume effects 

In order to investigate the presence of discretization effects and to study the renormaliza- 
tion scale dependence of the RCs of four-fermion operators we follow a procedure close to 
the one described in sects. 13. II and 13.31 for the case of bilinear quark operators. 
We start by constructing the RGI combinations 

Z^^^ = C-\^^)Z(^^) (34) 
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for both the PC and PV operators, where the evolution matrix C{fi), for the complete 
basis has been computed, in the RI/MOM scheme, at the NLO in perturbation theory [.45] . 

We then rescale the RGI combinations, computed at the different values of the lattice 
spacing, to a common reference scale a, which in this case we choose to be the value of 
the lattice spacing at /? = 6.45. To that purpose, we compute the renormalization scale 
independent ratios 

R{a, a) = Z^^\a)-^Z^^\a) = Z{a, ^i)-^Z{a, /i) . (35) 

which are the analogous of those defined in eq. Q for bilinear operators. We then use 
these ratios to rescale the RGI RCs at the three values of the couphng different from 
j3 = 6.45. 

The rescaled RGI combinations should be independent of both the renormalization 
scale and the lattice spacing. This is true, however, only up to discretization effects and 
higher-order perturbative corrections not taken into account in the NLO perturbative 
evaluation of the function C{n). A comparison of the rescaled RGI combinations obtained 
at the same value of the renormalization scale, but at different values of the lattice spacing, 
provides an estimates of discretization effects. On the other hand, by studying at fixed 
lattice spacing the dependence of the RGI combinations on the renormalization scale we 
can investigate the effect of higher orders perturbative contributions not included in the 
evaluation of the evolution matrix C(/i). 

In fig. ^1 we show the numerical results for the rescaled RGI combinations, in both 
the PV and PC function of the renormalization scale. In these plots we only 

show the results for the plus sector. For the minus sector the situation is similar. As can 
be seen from the plots, most of the matrix elements of both Z^'-'^ and Z^'^^ do not show 
any significant dependence on the renormalization scale. There are three cases, however, 
namely the matrix elements 33, 44 and 55, in which the plateau is worse than in the 
others. Since deviations from the expected constant behaviour are observed in a similar 
way at all values of the lattice spacing, and the results obtained at different f3 are in 
good agreement among each others, we conclude that lattice artifacts are quite small and 
interpret these deviations as mainly due to the effect of N^LO perturbative corrections not 
included in the evaluation of the function C(/i). This explanation is also supported by the 
fact that the operators Qs, Q4, and Q3, Q4, Q5 have very large anomalous dimensions 
at the leading order. We emphasize that a better control of the renormalization scale 
dependence of these operators could be achieved either by a perturbative calculation of 
the N^LO anomalous dimensions or, even better, by a non-perturbative study of the 
running by using an iterative matching technique which involves several lattice scales. 
This approach, which has been already implemented within the context of the SF scheme 
(see for instance ref. [12] for a study of the pseudoscalar RC Zp) can be implemented 
with other NPR techniques as well, like the RI/MOM or the x-space method. 

In fig.^Jwe show the results for the mixing coefficients A as a function of the renormal- 
ization scale /i. At fixed lattice spacing, these quantities are expected to be independent 
of fi. Indeed, in most of the cases we observe from the plots a reasonably good scale 
independence at large p"^, in particular for those coefficients of relatively large magnitude. 
The systematic error induced by the effect of a residual scale dependence is accounted for 
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Figure 10: Results for some of the RGI combinations as obtained at the four values of the 
lattice coupling and rescaled to j3 = 6.45. 
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Figure 11: Results for some of the mixing coefRcients Aij at the four values of lattice coupling. 



in the final determination of tlie RCs. 

Mixing among operators of different naive cliirality, parameterized by tlie matrix A, 
is a consequence of tlie explicit chiral symmetry breaking induced by the Wilson term. 
Therefore this mixing is expected to disappear in the continuum limit. To verify this 
expectation, we show in the left plot of fig. ^1 the results for some of the matrix elements 
of A as a function of the lattice spacing. We see from this plot that the absolute values 
of the As decrease as the continuum limit is approached. In the right panel of fig. ^1 we 
present the same matrix elements of A divided by Qpiaq, where gpi^^ = 5'^/(|Tr(f/p)) is 
the boosted coupling defined from the plaquette. If the perturbative expansion in terms 
of 

9piaq is rapidly convergent, the ratios ^/g^iaq are expected to be flatter than the As 
themselves. The plot in fig.ll2lshows that this is indeed the case. The residual dependence 
on the lattice spacing, observed in the figure, signals the presence either of higher orders 
in the perturbative expansion or of finite lattice spacing effects. 

To conclude the analysis of the systematic errors which may affect the determination 
of the RCs of four-fermion operators we investigate the presence of finite volume effects. 
As already done for bilinear quark operators (see sec. we compare the results for the 
RCs obtained at /3 = 6.0 on two different volumes, namely 16'^ x 52 and 24'^ x 64. The 
parameters of the simulation on the smallest volume are those given in table H For the 
calculation on the largest volume, we have used a set of 480 gauge field configurations and 
computed quark propagators at four values of the Wilson parameter, namely k =0.13180, 
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Figure 12: Results for some of the matrix elements of A (left) and for the ratios A/g^i^^ (right), 
where g^i^^ = / (^Tr(f/p)) is the boosted coupling defined from the plaquette, as a function of 
the lattice spacing. 



0.13280, 0.13770, 0.13440 (ncr = 0.135220(6)). The comparison of the results is illustrated 
in fig. ^Jfor some of the matrix elements of the RCs and mixing coefficients. We see from 
the plots that, in some cases, the presence of finite volume effects is visible, though they 
become negligible at large values of the renormalization scale [fi ^ 2.5 GeV). Since our 
final determination of the RCs is obtained from a fit in the interval fi = [2.5, 3.0] GeV (see 
below), we find that the results obtained on the different volumes are compatible within 
the statistical and the estimated systematic errors. For this reason, we do not include in 
the final error any additional uncertainty due to final volume effects, although this point 
would require more investigations. In order to better illustrate this comparison, when 
quoting our final results (see tables |31|7j) we will present at /3 = 6.0 the determinations 
obtained on both the lattices with different size. 



5.4 Results 

In order to obtain the RCs at a given scale we proceed in the following way. We extract 
the RGI combinations Z^^^ and Z^'^^ by fitting to a constant the plateau in fig. ^[ and 
similarly for the other matrix elements, in the interval fi = [2.5, 3.0] GeV. We then use 
the renormalization group evolution at the NLO to obtain the value of the RCs at the 
desired scale. A standard choice of the scale, at which matrix elements are matched with 
Wilson coefficients in the effective weak Hamiltonian, is /i = 2 GeV. Since this value lies 
in most of the cases within the range of momenta at which the RCs are directly computed, 
we estimate the systematic error in the following way: we take the RCs computed non- 
perturbatively at a given scale /io and run them up or down to the desired scale /x = 2 GeV. 
The systematic error is then estimated from the deviations of the results obtained by 
performing different choices of /iq. The same procedure has been also applied for the 
determinations of the mixing coefficients As. The only difference is that, being free of 
anomalous dimension, these quantities do not evolve with the renormalization scale. 

We emphasize that the systematic error induced by the renormalization group running 
of the RCs should be regarded as an uncertainty related mainly to the perturbative 
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Figure 13: Comparison for some of the matrix elements of the RCs and mixing coefficients as 
obtained on two different volumes at (3 = 6.0. 



estimates of the Wilson coefficients of the effective Hamiltonian, rather than to the non- 
perturbative the determination of the RCs. Indeed, had we chosen to quote the values of 
the RCs at the renormalization scale at which they are directly determined in the lattice 
calculation, this error would be absent. 

Our final estimates of the RCs of four-fermion operators in the RI/MOM scheme at 
the scale /i = 2 GeV are collected in tablesEJElElandlTlfor the PC^ and PV^ respectively. 
As mentioned before, for the case of /3 = 6.0, we present the results obtained on lattices 
with two different sizes. 



Conclusions 

In this paper we have presented the results of an extensive lattice calculation of the RCs 
of bilinear and four-quark AF = 2 and AI = 3/2 operators by using the RI/MOM 
NPR method. Several sources of systematic errors, including discretization errors and 
finite volume effects, have been investigated. When possible, we have also compared the 
results with those obtained by using other non-perturbative approaches, like the WI and 
the SF methods. In these cases, we find an agreement which is typically at the level 
of 1%. This comparison supports the conclusion that the RI/MOM method allows to 
obtain an accurate non-perturbative determination of the RCs of lattice operators. At 
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/5 


6.0 (16^ X 52) 


6.0 (243 X 64) 


6.2 


6.4 


6.45 


^11 


0.609(4)(13) 


0.604(3)(4) 


0.634(5)(5) 


0.671(10)(9) 


0.682(7)(5) 


Ai3 

Am 

Ai5 


-0.067( 1)( 4) 
-0.021( 1)( 1) 
0.010( 1)( 3) 
0.004( 1)( 3) 


-0.068( 2)( 4) 
-0.022( 1)( 1) 
0.014( 1)( 3) 
0.003( 1)( 2) 


-0.058( 2)( 4) 
-0.016( 1)( 2) 
0.011( 2)( 4) 
0.006( 1)( 1) 


-0.048( 1)( 3) 
-0.012( 1)( 2) 
0.013( 1)( 3) 
0.004( 2)( 1) 


-0.045( 2)( 2) 
-0.012( 1)( 2) 
0.013( 1)( 3) 
0.006( 1)( 3) 


■^22 
■^23 


0.669(3)(11) 
0.274(2)(8) 


0.662(3)(3) 
0.276(2)(2) 


0.691(3)(5) 
0.298(4)(7) 


0.732(8)(10) 
0.323(8)(9) 


0.740(8) (7) 
0.328(6)(10) 


A21 
A24 
A25 


-0.050( 1)( 6) 
-0.199( 4)( 7) 
0.020( 1)( 2) 


-0.053( 1)( 5) 
-0.194( 3)( 3) 
0.019( 1)( 1) 


-0.047( 2)( 5) 
-0.158( 6)( 7) 
0.013( 1)( 3) 


-0.040( 1)( 4) 
-0.111( 2)(12) 
0.008( 1)( 3) 


-0.038( 1)( 3) 
-0.102( 4)(17) 
0.008( 1)( 4) 


^32 
^33 


0.029(1)(1) 
0.390(3)(31) 


0.028(1)(1) 
0.376(2)(15) 


0.025(1)(1) 
0.357(4)(19) 


0.024(2)(1) 
0.344(5)(12) 


0.022(1)(1) 

0.331(6)(18) 


A31 
A34 
A35 


0.014( 1)( 4) 
0.269( 6)(18) 
-0.009( 0)( 1) 


0.015( 1)( 3) 
0.260( 4)( 3) 
-0.008( 1)( 1) 


0.015( 1)( 3) 
0.215( 7)( 7) 
-0.006( 1)( 1) 


0.013( 1)( 1) 
0.153( 3)(16) 
-0.003( 1)( 2) 


0.013( 0)( 3) 
0.131( 7)(24) 
-0.002( 1)( 5) 


Z44 
^45 


0.478(3)(33) 
-0.024(1)(1) 


0.463(2)(15) 
-0.023(1)(1) 


0.445(4) (20) 
-0.019(1)(1) 


0.434(5)(14) 
-0.017(2)(2) 


0.419(5)(21) 
-0.015(2)(1) 


A41 
A42 
A43 


0.005( 0)( 2) 
0.007( 0)( 2) 
0.174( 4)( 6) 


0.005( 0)( 2) 
0.006( 0)( 0) 
0.164( 2)( 3) 


0.006( 0)( 2) 
0.006( 1)( 0) 
0.136( 4)( 7) 


0.005( 0)( 1) 
0.005( 1)( 2) 
0.100( 2)(10) 


0.005( 0)( 1) 
0.003( 1)( 3) 
0.085( 5)(11) 




-0.253(2)(13) 
0.786(4)(15) 


-0.257(2)(6) 
0.788(3)(9) 


-0.291(3)(10) 
0.852(5)(19) 


-0.329(6)(13) 
0.934(12)(23) 


-0.337(6)(13) 
0.953(9)(19) 


A51 
A52 
A53 


0.005( 0)( 3) 
0.013( 0)( 1) 
0.107( 2)( 2) 


0.005( 1)( 2) 
0.012( 0)( 1) 
0.100( 2)( 3) 


0.007( 1)( 2) 
0.009( 1)( 1) 
0.081( 2)( 6) 


0.007( 1)( 1) 
0.006( 1)( 2) 
0.060( 1)( 6) 


0.007( 0)( 2) 
0.006( 1)( 2) 
0.052( 3)( 6) 



Table 4: Results for the RCs of the PC+ sector in the RI/MOM scheme at the scale ijl=2 GeV. 
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-0.034( 1)( 5) 
0.007( 1)( 1) 
-0.013( 1)( 6) 
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-^22 
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-0.277(2)(10) 


0.669(2)(3) 
-0.282(2)(5) 


0.700(4)(5) 
-0.306(4)(9) 


0.743(9)(11) 
-0.335(8)(12) 


0.749(8)(7) 
-0.340(5)(10) 


A21 
A24 
A25 


-0.045( 1)( 3) 
0.138( 2)( 6) 
0.012( 1)( 3) 


-0.050( 2)( 3) 
0.134( 3)( 1) 
0.012( 1)( 1) 


-0.040( 2)( 3) 
0.113( 3)( 5) 
0.010( 1)( 2) 


-0.033( 1)( 1) 
0.078( 2)( 8) 
0.007( 1)( 4) 


-0.028( 2)( 4) 
0.064( 5)(13) 
0.005( 1)( 5) 


^32 
^33 


-0.027(1)(1) 

0.402(3)(30) 


-0.028(1)(1) 

0.386(2)(14) 


-0.023(1)(1) 
0.366(4)(19) 


-0.021(2)(1) 
0.349(5)(13) 


-0.021(1)(2) 
0.335(7)(19) 


A31 
A34 
A35 


-0.008( 1)( 5) 
0.203( 3)( 9) 
0.000( 0)( 6) 


-0.011( 1)( 1) 
0.195( 3)( 3) 
0.000( 1)( 2) 


-0.007( 1)( 2) 
0.163( 4)( 9) 
-0.002( 1)( 2) 


-0.003( 1)( 3) 
0.115( 3)(13) 
-0.002( 1)( 4) 


-0.004( 2)( 5) 
0.094( 8)(18) 
-0.005( 2)( 9) 


Z44 
■^45 


0.309(5)(44) 
-0.007(1)(3) 


0.300(3)(15) 
-0.007(1)(1) 


0.280(5)(20) 
-0.016(1)(2) 


0.258(5)(16) 
-0.020(2)(1) 


0.238(7)(37) 
-0.019(4)(18) 


A41 
A42 
A43 


0.004( 2)(32) 
-0.010( 1)( 7) 
0.331(18)(85) 


0.000( 1)( 1) 
-0.012( 1)( 2) 
0.303( 5)( 6) 


0.002( 1)( 4) 
-0.004( 1)( 4) 
0.248( 9)( 5) 


0.003( 2)( 4) 
0.000( 2)(13) 
0.174( 4)(19) 


-0.021(14)(63) 
-0.001(10)(55) 
0.165(22)(184) 




0.120(2)(3) 
0.713(3)(10) 


0.124(2)(1) 
0.712(3)(5) 


0.133(3)(3) 
0.745(4)(16) 


0.139(5)(2) 
0.800(12)(18) 


0.133(7)(15) 
0.804(13)(26) 


A51 
A52 
A53 


-0.003( 1)( 6) 
0.013( 1)( 1) 
-0.113( 3)( 5) 


-0.005( 1)( 4) 
0.013( 1)( 1) 
-0.107( 2)( 2) 


-0.010( 1)( 3) 
0.009( 1)( 2) 
-0.087( 3)( 5) 


-0.014( 1)( 4) 
0.007( 1)( 2) 
-0.064( 1)( 6) 


-0.014( 2)( 4) 
0.006( 1)( 4) 
-0.058( 2)( 6) 



Table 5: Results for the RCs of the PC sector in the RI/MOM scheme at the scale fi=2 GeV. 
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Zu 
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0.425(5)(16) 


0.417(5)(22) 


Z45 


-0.024(1)(1) 


-0.022(1)(1) 


-0.019(1)(1) 


-0.016(2)(1) 


-0.015(2)(3) 


■^54 


-0.255(2)(11) 


-0.260(2)(5) 


-0.295(3)(12) 


-0.333(6)(14) 


-0.339(6)(13) 


Z55 


0.787(4) (15) 


0.790(3) (9) 


0.854(5)(19) 


0.936(12) (24) 


0.952(10)(19) 



Table 6: Results for the RCs of the Py+ sector in the RI/MOM scheme at the scale ii=2 GeV. 
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-0.007(1)(1) 
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0.126(1)(4) 
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0.135(3)(3) 


0.140(5)(1) 


0.130(10)(36) 


-^55 


0.715(3)(12) 


0.714(3) (7) 


0.746(4) (17) 


0.802(12)(18) 


0.806(16) (32) 



Table 7: Results for the RCs of the PV sector in the RI/MOM scheme at the scale fi=2 GeV. 
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the same time, the method is extremely simple to implement and can be applied to a 
large class of operators. The only notable exception is represented by those operators the 
renormalization of which requires subtraction of power divergences. In the latter case, 
and in particular in the case of the AI = 1/2 four-fermion operators, the use of non 
gauge invariant correlation functions renders the applicability of the RI/MOM method 
extremely difficult, if not impossible in practice. In order to renormalize these operators 
the gauge invariant NPR method in x-space [I]-|nj might be the method of choice. 

Acknowledgments 

It is a pleasure to thank C.J.D. Lin, A. Le Yaouanc, CT. Sachrajda and M. Testa for 
many interesting discussions on the subject of this paper. 



Appendix A 

In this appendix we show that the RCs of bilinear quark operators obtained with the 
RI/MOM method in the chiral limit, at sufficiently large values of the external momentum 
and at zero momentum transfer, are automatically improved at 0{a). 

The renormalized and improved version of the bilinear operator Or = qTq has the 
form jlH]: 

d[ = ZrOl = ZriqTq + acr04,r + ac'rEr) (36) 

where Zr is the multiplicative RC in which we also include here the 0{a) linear dependence 
on the quark mass, Z-p = Zp (1 + fepam). The coefficients b^, cp, c'r express the mixing of 
Or with higher-dimension operators. This mixing has to be taken into account in order 
to improve the operator at 0{a). The operator Er in eq. (j36|l is given by: 

Er = q [t(^ + mo) + + mo)r] q (37) 

where p + mo is a shorthand for the entire lattice fermion operator, including the Wilson 
and the SW-clover terms. Therefore E-p vanishes by the equation of motion and it only 
contributes to contact terms when inserted in correlation functions. The operator 04,r 
in eq. ()36|) is a gauge-invariant, dimension-four operator which does not vanish by the 
equation of motion. In the cases of the vector, axial-vector and tensor operators, 04,r is 
given by 

Oiy = (qa^^^q) 

0,,A = d.iq^'q) (38) 
04,T = df, (q-f^q) - du {q-ff,q) 

while this mixing is absent for the scalar and pseudoscalar densities. 

In order to improve the off-shell correlation functions, we also consider the renormal- 
ized improved quark and antiquark fields ||46], 

q = Zg^^"^ 1 + acq{p + mo) + acj^Gi^ q (39) 
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l + ac'{-p +mo) -acNGi' 



with Zq = Z^l (1 + hqam). In terms of these fields we can define the renormahzed improved 
quark propagator: 

S{p) = J d^xe-'P-'' (g(x)f(O)) , (40) 
which, by using eq. (jSHl), can be related to the corresponding lattice quantity S{p) by 

S{p) = [S{p) + 2a c'q + 2a cngi ^^ S{p)] (41) 

up to C(a^) terms. 

Let us now discuss the effect of the off-shell improvement, defined by eqs. (jHUj) and (jSHl), 
in the determination of the RCs with the RI/MOM method. The relevant renormahzed 
improved correlation function is defined in terms of improved operators and external fields: 



I d'^x d% e-'P-^^'P'-y (g(x) {qTq + a Cr0,,r + a c'rEr), f (y)) . (42) 

The RI/MOM method consists in imposing that the forward amputated Green function, 

Ai{p)=S{p)-'Gi{p,p)S{p)-\ (43) 

computed in a fixed gauge and renormahzed at a given scale = fi"^, is equal to its 
tree-level value. In practice, this condition is implemented by requiring: 

Zr Tiip) = Zr Tr [Ai{p) Pr] \p2=,^ = 1 (44) 

where Pr is the Dirac projector. 

The Green functions Gf (p) and Af (p) can be easily related to their unimproved coun- 
terparts. By using eqs. (jHUjl - ljl^ . and up to terms of O(a^), one finds: 

G'Ap) = Grip) + a cr G,,r{p) + a (c'r + c',) (r S{p) + S{p)t) + 

0- Cngi UGlip) + G'rip)^i>) (45) 



where Gr{p) is the Green function of the unimproved operator Or, constructed with 
unimproved external quark fields, and G'4,r(p) is the analogous quantity for the operator 
O4J defined in eq. fl36p . The improved function Af (p) is then obtained from Gf(p) by 
amputating the external legs with improved quark propagators, according to eq. (021) • 
The result is: 

A^;(p)=Ar(p) + acrA4,r(p) + a{c'r + c'q) (t S{p)-' + S{p)-'r) - 

2ac'qZ-' {S{p)-'A'Ap) + A'Ap)S{pr') - ac^Gl{^Mkp) + AUp)^^) (46) 
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Finally, by projecting eq. pUj) onto the tree-level form factor with the projector Pp? one 
obtains: 

riip) = Trip) + a cr T,^r{p) + 2 a (c'r + c',)Tr [S{p)-'] - 
2a 4 Z-^Tr [S{p)-' Ai{p) Pr - A'Ap) S{p)-' Pr) + 
a CNGiTr {i^ {p) Pr + A^^ {p) i^Pr) . (47) 

The various contributions to the amputated projected Green function in eq. (jlTj) are easily 
identified. The third term on the r.h.s. of eq. ()47p comes from the operators vanishing by 
the equation of motion and it is proportional to the trace of the improved inverse quark 
propagator. Since at large values of p^ this quantity is proportional to the quark mass, 
such a term does not affect the determination of the renormalization constants in the 
chiral limit. 

One may be easily convinced that also the last two terms in eq. (14 7|) proportional to c'q 
and cngi respectively, vanish in this limit. We first observe that, for vanishing quark mass, 
the improved propagator S{p)~^ at large p^ is proportional to Therefore, the two terms 
have the same form. In addition, we note that contributions from Af (p) proportional 
to the tree-level form factor, F, vanish because they reduce to Tr(j^). For the same 
reason, one also finds the vanishing of all possible contributions coming from different 
(dimensionless) form factor which do not depend on the quark mass. For instance, in the 
case of the vector current, the Green function A^(p) also contains a term proportional 
to p^ii/p^i which again gives a contribution vanishing as Tr (j^) to the last two terms of 
eq. fj47|) . A non vanishing contribution may come from a form factor proportional to 
m|^7^/p^, but then this contribution vanishes in the chiral limit. 

In the case of the vector, axial-vector and tensor operators, on-shell (9(a)-improvement 
also requires to consider the mixing of the bare operators Or with the dimension-4 oper- 
ators 04^r5 whose contribution is represented by the second term on the r.h.s. of eq. (jTTj) . 
For the scalar and pseudoscalar densities this mixing is absent. In the other cases, the 
operators 04,r are listed in eq. (jHHj) . A simple observation is that the operators O^^y all 
have the form of a 4-divergence. Therefore, their forward correlation functions, r4^r(p) in 
eq. (I47j) . vanish identically. This completes the proof that only the unimproved correlator 
Ty{p) survives at large p^ and in the chiral limit on the r.h.s. of eq. (j47p and contributes to 
the calculation of the RCs with the RI/MOM method. Such a calculation can thus be per- 
formed even in the lacking of a non-perturbative determination of the coefficients cr and 
c'r. Note also that, although the correlation function rf(p) is obtained from G^{p,p) by 
amputating the external legs with improved quark propagators, in practice, in the chiral 
limit, improving the quark propagator is also an unnecessary step for the determination 
of the RCs. 

Appendix B 

In this appendix we show that, in the case of four-fermion operators, both single and 
double Goldstone boson poles can appear in the chiral limit in the relevant correlation 
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functions, and may thus affect the RI/MOM NPR procedure at finite values of the external 
momenta. For convenience, we work with non vanishing quark masses and we will consider 
the zero quark mass limit only at the end of the calculation. 

Let us consider the Green function of a four fermion operator Oi with four external 
quark states and two different external momenta 

Gi{p, p') = J d^xi d^X2 d^X3 d\i e-^^'(^^+^«)+*P'(^2+^*) X 

= {Q\T[^-{p)^\p')Oi{Q)^-{p)^'^{p')m (48) 
By applying the LSZ reduction formula we can write 



G,{p,p') = mrwnpViq))- — ^{7r'''{q)\T[o,io)r{p)np')m + 



NP 



{0\T[i>-{p)il;b{p')O,i0)^-ip)i>'iip')]\0) + ... (49) 



where q^ = p^ — p'^ and we have identified the Goldstone boson with the vr meson. In 
eq. ()49|) the symbol NP[. . .] indicates the part of the correlation function without pions in 
the intermediate states and the dots represent terms with pion propagators corresponding 
to pions different from |7r''"(g)) (i.e. |vr'^"(q')), \T^'^'^{q)) and |7r'"^(g))). If we now apply the 
reduction formula again we obtain 

G,{p,p') = {0\T[rm\p')]\n''-iq))-^^ X 



(7r^"(g)|O.(0)|7r'^'^(g))^ ^(7r^'^(g)|T[^^(0)^^(p')]|0) + 

'-'ii ^ 



q^ — m. 



7c'^{q)\T[omr{pmp')m]} + . . . . m 

In the implementation of the RI/MOM method we have considered the case = p'^, 
corresponding to q^ = 0. From eq. (j5(jp it is then clear that both single and double 
pion poles can appear in the amputated Green functions in the chiral limit. In the case 
of parity violating operators, however, since by parity (7r^"(0)|Of^(0)|7r^'^(0)) = 0, only 
single poles can be present. 
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